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Abstract 

In this paper we investigate the saturation physics in diffractive deep inelastic electron-ion scat- 
tering. We estimate the energy and nuclear dependence of the ratio a dt ff ja tot and predict the 
xjp and [3 behavior of the nuclear diffractive structure function F^^\Q 2 , (3 , x p) . Moreover, we 
analyze the ratio R^f^iQ 2 ^ ^ x p) = F2A1 I F2A2 > wri ich probes the nuclear dependence of the 
structure of the Pomeron. We show that saturation physics predicts that approximately 37 % of 
the events observed at eRHIC should be diffractive. 
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I. INTRODUCTION 



One of the frontiers of QCD which are intensely investigated in high energy experiments is 
the high energy (small x) regime, where we expect to observe the non-linear behavior of the 
theory. In this regime, the growth of the parton distribution should saturate, forming a Color 
Glass Condensate (CGC). (For recent reviews see, e.g. Refs. 

mm 

. In fact, signals of 

parton saturation have already been observed both in ep deep inelastic scattering at HERA 
and in deuteron-gold collisions at RHIC (See, e.g. Ref. [5]). However, the observation of this 
new regime still needs confirmation and so there is an active search for new experimental 
signatures. Among them, the observables measured in diffractive deep inelastic scattering 
(DDIS) deserve special attention. As shown in Ref. Q], the total diffractive cross section is 
much more sensitive to large-size dipoles than the inclusive one. Saturation effects screen 
large-size dipole (soft) contributions, so that a fairly large fraction of the cross section is 
hard and hence eligible for a perturbative treatment. Therefore, the study of diffractive 
processes becomes fundamental in order to constrain the QCD dynamics at high energies. 

Significant progress in understanding diffraction has been made at the ep collider HERA 
(See, e.g. Refs. P, ls[ El). Currently, there exist many attempts to describe the diffractive 
part of the deep inelastic cross section within pQCD (See, e.g. Refs. 

QQ)- One of 

he most successful approaches is the saturation one :6| based on the dipole picture of DIS 







14j . It naturally incorporates the description of both inclusive and diffractive events 



in a common theoretical framework, as the same dipole scattering amplitude enters in the 
formulation of the inclusive and diffractive cross sections. In the studies of saturation effects 
in DDIS, non-linear evolution equations for the dipole scatte ring amplitude have been derived 

and the charm contribution 
Q,~" t 7a t ~e not y e t precise enough.no, 



00, EL 

Jm|, new measurements proposed 



estimated 23^. However, as shown in Ref. 
do they extend to sufficiently small values of xjp, to discriminate between different theoretical 
approaches. 

Other source of information on QCD dynamics at high parton density is due to nuclei 



which provide high density at comparatively lower energies. Recently, in Ref. 



2J, we 



have estimated a set of inclusive observables which could be analyzed in a future electron- 



ion collider 



Our results have demonstrated that the saturation physics cannot be 



disregarded in the kinematical range of eRHIC. Our goal in this work is to understand 
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to what extend the saturation regime of QCD manifests itself in diffractive deep inelastic 
eA scattering. In particular, we will study the energy and nuclear dependence of the ratio 
between diffractive and total cross sections (o'diff/o'tat)- HERA has observed that the energy 
dependence of this ratio is almost constant for different mass intervals of the diffractively 
produced hadrons over a wide range of photon virtualities Q 2 j^j]. This ratio is to a good 
approximation constant as a function of the Bjorken x variable and Q 2 . Moreover, we make 



predictions for more detailed diffractive properties, such as t 



rose embodied in the diffractive 



structure function F 2 D ^ (Q 2 , (3,Xjp). Motivated by Refs. 

of the ratio between nuclear diffractive structure functions R d 2(f A2 (Q 2 , (3, x P ) = F 2 D $/E 
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271 ] we also analyze the behavior 



D(3) / p D(3) 
2,A2 > 



where A\ and A 2 denote the atomic number of the two nuclei. It is important to em pha size 



that diffractive processes in eA collisions were studied in Refs. 
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32, 




351 ] . Here we extend these studies for a large number of observables, considering the dipole 



approach and a generalization for nuclear targets of the CGC dipole cross section proposed 
in Ref. [sf|. As this model successfully describes the HERA data, we believe that it is 
possible to obtain realistic predictions for the kinematical range of the electron-ion collider 
eRHIC. 

This paper is organized as follows. In next Section we present a brief review of the 
dipole picture. We present the main formulae for the dipole cross section and the diffractive 
structure function. In Section IIHI we introduce the overlap function for diffractive events, 
which allows us to find out the average dipole size which contributes the most to this process, 
and analyze its dipole size and nuclear dependences. Moreover, we estimate the different 
contributions to the diffractive structure function and present our predictions for F 2 ^ and 
R2i A2- Finally, in Section HVl we summarize our main conclusions. 



II. DIPOLE PICTURE OF DIFFRACTIVE DIS 

In deep inelastic scattering, a photon of virtuality Q 2 collides with a target. In an appro- 
priate frame, called the dipole frame, the virtual photon undergoes the hadronic interaction 
via a fluctuation into a dipole. The wavefunctions \ipr\ 2 and \4>l\ 2 , describing the splitting 
of the photon on the dipole, are given by Q] : 

\M^r)\ 2 = 3 -^J2e 2 4Q 2 a 2 (l-a) 2 K 2 (er) (1) 
n f 
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hM«,r)| 2 = ^^e){[a^{l-afVKl{er)+m)Kl{er)} (2) 

for a longitudinally and transversely polarized photon, respectively. In the above expressions 
e 2 = a(l — a)Q 2 +m 2 , i^" an d i^i are modified Bessel functions and the sum is over quarks of 
flavor / with a corresponding quark mass m/. As usual a stands for the longitudinal photon 
momentum fraction carried by the quark and 1 — a is the longitudinal photon momentum 
fraction of the antiquark. The dipole then interacts with the target and one has the following 
factorized formula [13] 

ct2*t(x,Q 2 ) = J dad 2 r \4>L,T(a,r)\ 2 a di p(x,r) . (3) 



Simil arly the total diffractive cross sections take on the following form (See e.g. Refs. 
1) 



bimilarly 

la 3,13; 

J — oo Cut j_q J 

where 



da D 

(4) 

t=o 



T\L_ 

B D dt 



d °T,L 



t=0 16 



l — J d 2 v jf tZa|* r ,L(a,r)| 2 (7jp(a;,r 2 ) , (5) 



dt 

and we have assumed a factorizable dependence on t with diffractive slope Bp. 

The diffractive process can be analyzed in more detail studying the behavior of the 
diffractive structure function F® {Q 2 , (3,xjp). In Refs. DO the authors have derived 
expressions for directly in the transverse momentum space and then transformed to 

impact parameter space where the dipole approach can be applied. Following Ref. [6( we 
assume that the diffractive structure function is given by 

F 2 D(3) (Q 2 , 0, xjp) = F° L + F° T + F° gtT (6) 

where T and L refer to the polarization of the virtual photon. For the qqg contribution 
only the transverse polarization is considered, since the longitudinal counterpart has no 
eadin g lo garithm in Q 2 . The computation of the different contributions was made in Refs. 



l37| and here we quote only the final results: 



30 6 r 1 / 2 

XlP F g l L (Q 2 , [3, xjp) = 32n , pBD E e}2 <W(1 - af%, (7) 

30 4 r 1 / 2 r i 

x F F° T (Q 2 ,P,x F ) = U8 ^ (3Bd E4 2 J aQ daa(l-a){e 2 [a 2 + (l-a) 2 ]^ 1 + m 2 ^ } 

(8) 



where the lower limit of the integral over a is given by «o 
introduced the auxiliary functions 



4mJ 



0,1 



rdrK 0:1 {er)a d ip{x]p, r)J 0t i{kr] 



and we have 



(9) 



For the qqg contribution we have 

81(3a s 



(l-z)Q 2 



dk 2 \n 



' 512tt 5 B d 
k 2 t ) 
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dz 



udu a dip (u/k t ,X]p)K2 



+ - 



(10) 



-U 



Ji(u) 



We use the standard notation for the variables (3 = Q 2 /(M X + Q 2 ), x F = (M| + Q 2 )/{W 2 + 
Q 2 ) and x = Q 2 /(W 2 + Q 2 ) = fixjp, where Mx is the invariant mass of the diffractive 
system and W the total energy of the -y*p (or 7*A ). When extending (J7J), (JHJ) and (jTUjl to 
the nuclear case we need to change the slope to the nuclear slope parameter, Ba- In what 

R 2 

follows we will assume that Ba may be approximated by Ba = -f-, where Ra is given by 
R A = 1.2 A 1 / 3 fm 



39] 



In the dipole picture the behavior of the total inclusive and diffractive cross sections, as 
well as the diffractive structure functions, is strongly dependent on the dipole cross section, 
which is determined by the QCD dynamics. Consequently, in the dipole picture the inclusion 
of saturation physics is quite transparent and straightforward, as the dipole cross section is 
closely related to the solution of the QCD non-linear evolution equations (For recent reviews 
see, e.g. Refs. QIISQ) 



a dip (x,r) = 2 J d 2 bAf(x,r,b) 



(11) 



where M is the quark dipole-target forward scattering amplitude for a given impact pa- 
rameter b which encodes all the information about the hadronic scattering, and thus about 
the non-linear and quantum effects in the hadron wave function. In what follows we will 
disregard the impact parameter dep endence [a dip = (JqN{x, r)\ and consider the phenomeno- 
logical model proposed in Ref. 36], in which a parameterization of J\f(x, r) was constructed 
so as to reproduce two limits analytically under control: the solution of the BFKL equa- 
tion for small dipole sizes, r 1/Q s (x), and the Levin- Tuchin law for larger ones, 
r ^> 1/Q s (x). Here, Q s denotes the saturation momentum scale, which is the basic quantity 
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characterizing the saturation effects, being related to a critical transverse size for the uni- 
tarization of the cross section, and is an increasing function of the energy [Q 2 = Qq (^t) a ]. 
A fit to the structure function ^(x, Q 2 ) was performed in the kinematical range of interest, 
showing that it is not very sensitive to the details of the interpolation. The dipole-target 
forward scattering amplitude was parameterized as follows, 

J ln(2/VQ 3 ) 



M(,,r) = {^m V ' ■*» •<?.<*)< 2, (12) 

1 - exp -*^ 2 (brQ s ) ; for r Q s (x) > 2, 

where a and b are determined by continuity conditions at rQ s (x) = 2, 7 S = 0.63, k = 9.9, 
A = 0.253, Ql = 1.0 GeV 2 , x = 0.267 x 10" 4 and Af = 0.7. Hereafter, we label the model 
above by IIM. The first line from Eq. (|12jl describes the linear regime whereas the second 
one describes saturation effects. When estimating the relative importance of saturation we 
will switch off the second line of (J 12)) and use only the first. This is a relevant check, since 
some observables may turn out to be completely insensitive to non-linear effects. Following 



Ref. |2j], we generalize the IIM model for nuclear collisions assuming the following basic 
transformations: a — > (Tq = A% x <7o and Q 2 s (x) ^ Q 2 A = x Q 2 (x). As already empha- 
sized in that reference, more sophisticated generalizations for the nuclear case are possible. 
However, as our goal is to obtain a first estimate of the saturation effects in these processes, 
our choice was to consider a simplified model which introduces a minimal set of assump- 
tions. In a full calculation we must use the solution of the BK equation, obtained without 
disregarding the impact parameter dependence as well as an initial condition constrained by 
current experimental data on lepton-nucleus DIS. 

III. RESULTS 

Before presenting our results for -fff- and for the diffractive structure function, let's 
investigate the mean dipole size dominating the diffractive cross section through the analysis 
of the photon-nucleus diffractive overlap function defined by 

H D (r,x,Q 2 ) = 2nr £ f da\%(a, r,Q 2 )\ 2 a% p (x,r,A). (13) 

i=T,L J 

In Fig. ^ we present the r-dependence of the photon-nucleus diffractive overlap function 
(normalized by A 2 ) at different values of the atomic number and Q 2 = 1 GeV 2 . A similar 



6 



3 











I ' 


1 1 1 1 1 1 






A = 2 (lull) 








A - 2 (linear) 








A= 197 (full) 








A = 197 (linear) 












_ 


Q 2 =1 








x = 10" 


i *■ 






IIM 























0.2 0.4 0.6 0.8 1 

r(fm) 



A = 2(full) 

A - 2 (linear) 

A ^ 197 (full) 



A = 197 (linear) 




FIG. 1: The r-dependence of the photon- nucleus diffractive overlap function (normalized by A 2 ) at 
different values of the atomic number and distinct saturation models: IIM (left panel) and GBW 
(right panel). 

analysis can be made for other values of Q 2 . The main difference is that at large values 
of Q 2 , the overlap function peaks at smaller values of the pair separation. For comparison 
he predictions obtained with a generalization of the Golec-Biernat Wusthoff (GBW) model 
f| for nuclear targets is also presented. When the saturation effects are included the IIM 
and GBW overlap functions present a similar behavior, strongly reducing the contribution 
of large pair separations. At large A only small pair separations contribute to the diffractive 
cross section. However, in the linear case, these two models present a very distinct behavior, 
which is directly associated with the different prescription for the linear regime. While 
the GBW model assumes that <7di V oc o§r 2 Q 2 s in this regime, the IIM one assumes o"^ p oc 
a [r 2 Q 2 ]' leff , where j e ff = 7 S + \n(2/rQ s )/ k XY is smaller than one. Firstly, it implies a 
different A dependence, since in the GBW model the product [croQ 2 ] 2 is proportional to A 2 
which cancels with the normalization term. In the IIM model we have [<ToQs 7e// ] 2 , which 

4 + 2 7 e // 

implies an i 3 dependence. When combined with the factor A which comes from the 

2(7 e //~l) 

normalization, we expect an A 3 dependence for the IIM overlap function. As 7^ < 1, 
the overlap function decreases for large A also in the linear regime. This behavior is seen in 
Fig. Secondly, in contrast to the GBW model which predicts a r 2 behavior for the dipole 
cross section in the linear regime, the IIM one leads to a r 2le ^ dependence. This different 
prescription for the r-dependence, when combined with the pair separation dependence of 
the wave functions, implies a strong modification on the contribution of the large dipole 
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FIG. 2: The ratio between the diffractive and total cross section as a function of x and W for 
different values of A and Q 2 . The black disk limit, &diff/otot = 1/2, is also presented. 

sizes, as observed in Fig. ^ It is important to emphasize that this contribution dominates 
the cross section, i.e. if we disregard the saturation effects, the diffractive cross section is 
dominated by soft physics. 

We now present a qualitative analysis of the A and energy dependence of the ratio 
(Jdiff/utot using the IIM model generalized for nuclear targets. Following Ref. p and 
assuming that a dip in the saturation regime can be approximated by ctq, the transverse part 
of the inclusive and diffractive cross sections, in the kinematical range where Q 2 > Q 2 , can 
be expressed as 

,4/Q* dr 2 r r 2 Q 2 dr 2 ( 1 \ r r 2 Q2 f00 dr 2 / 1 \ 

(14) 
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and 
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4/Q 2 ^2 2^2 



dr 2 



Q 2 r 2 



r 2 2 



dr 2 



4/Q1 



In order to obtain an approximated expression for the ratio we will disregard the in- 
dependence of the effective anomalous dimension, i.e. 7 e // = 7 = cte. In this case, we 
obtain adiff/ctot ~ [gt] 1-7 - Assuming 7 = 0.84, as in Ref. |36|], we predict that the ratio 
decreases with the photon virtuality and presents a weak energy dependence. However, an- 
alyzing the A- dependence, we expect a growth of approximately 30 % when we increase A 
from 2 to 197. In the kinematical range where Q 2 < Q 2 the ratio of cross sections presents a 
similar behavior. The main difference is that in the asymptotic regime of very large energies 
the cross section for diffraction reaches the black disk limit of 50% of the total cross section. 

In Fig. |2| we show the ratio (Jdiff/^tati computed with the help of © and (jlj), as a function 
of W and x for different values of A. The black disk limit, Udiff/^tot = 1/2, is also presented 
in the figure. We can see that the ratio depends weakly on W and on x but is strongly 
suppressed for increasing Q 2 . This suggests that in the deep perturbative re gion , diffraction 
is more suppressed. This same behavior was observed in diffractive ep data [26] . Moreover, 
the energy dependence of the ratio is remarkably flat, increasing with A, becoming 37 % (30 
%) larger for gold in comparison to proton (deuteron). This behavior agrees qualitatively 
with the previous calculation of jj^l and with our previous estimate. Similar results have 
been obtained in the pioneering work of Ref. j3| in a different context. The appearance of 
a large rapidity gap in 37 % of all eA scattering events would be a striking confirmation of 
the saturation picture. 

In Fig. IH1 we show our predictions for the diffractive structure functions 
xjpF^ 3 \x]p, (3,Q 2 ) as a function of (3 and different nuclei. We also present the linear pre- 
diction for X]pF^ 3 \ It is important to emphasize that a linear ansatz for the dipole cross 
section would not describe the HERA data. However, in order to estimate the importance of 
the saturation physics and clarify its contribution at different kinematical ranges, a compar- 
ison between these two predictions is valid. We can see that the normalization of xjpF®^ 1 
is strongly reduced increasing the atomic number, what is expected from our analysis of 
the diffractive overlap function. Moreover, although the photon wavefunction determines 
the general structure of the /9-spectrum p, |37(, the qqg component, which dominates the 
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FIG. 3: Diffractive structure function F 2 as a function of (5 and distinct nuclei. The qqg 
component of the diffractive structure function is explicitly presented. 




FIG. 4: Diffractive structure function F 2 as a function of and distinct nuclei. The transverse 
and qqg components of the diffractive structure function are explicitly presented. 

region of small (3, has its behavior modified by saturation effects and changes the behavior of 
xjpF 2 D ^ in this region. Moreover, the diffractive structure function becomes almost flat at 
intermediate values of (3 and large A. In Fig. Hjwe show an amplification of the lower curves 
in Fig. Eland include also the qqx component. In doing this another interesting feature of 
diffraction off nuclear targets emerges, namely, the relative reduction of the qqg component 
with respect to the qq one. 

In Fig. Elwe show our predictions for X]pF 2 ' (xp, (3, Q 2 ) as a function of xjp and different 
values of (3, Q 2 and A. Our choice for the combination of values of f3 and Q 2 was motivated 
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FIG. 5: Predictions for the diffractive structure functions xpF 2 (xp, (3, Q 2 ) plotted as a function 
of xjp for different values of 0, Q 2 and A. 



by the HERA results |26(. The xjp dependence comes from the dipole cross section, which 
in our case is given by the IIM model generalized to nuclear targets. We find that xpF 2 ^ 
increases at small values of xp. However, as the saturation scale grows with A, the xp 
becomes weaker when we increase the atomic number. 
In Fig. El we show our predictions for the ratio 



Ra{?A2(P>Q 2 i x 1p) 

as a function of (3 and xp. In our calculation we assume that A 2 = 2. Our analysis is 



F 2 y i \(3,Q 2 ,xp) 



F 2 D $(P,Q 2 ,xp) 
i we assume that 

motivated by Rets. j20. 27|. In these papers it was suggested that the nuclear dependence of 



auu x ip. 



this ratio can help us to establish the universality of the Pomeron structure. In particular, 
in Ref. 27], the authors have pointed out that if R A t ( A2 (j3,Q 2 ,xp) = 1 one can conclude 
that the structure of the Pomeron is universal and the Pomeron flux is A independent. On 
the other hand, if R^(f A2 (/3,Q 2 ,xp) = f(A 1 ,A 2 ), the structure is universal but the flux 
is A dependent. From our previous analysis we can anticipate that in the dipole picture, 
assuming the presence of saturation effects, this ratio will be A dependent. Moreover, its 
behavior will be determined by the saturation scale. In Fig. El we observe a strong decrease 
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FIG. 6: The ratio R d M f A2 (Q 2 , P,x P 
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of R2i A2 as a function of A. At the same time this ratio is remarkably flat at all values of A, 
xp and (3 > 0.2. However, at small (3, it presents a steeper dependence directly associated 
with the nuclear dependence of the qqg component of the diffractive structure function (See 
Fig. EJ). In order to estimate how much of the flat behavior is due to saturation we calculate 
Ra(a2 again using only the linear part of the dipole cross section, as discussed above, for the 
heaviest target (A = 197), for which the saturation effects are expected to be dominant and 
show our results in Fig. As it can be seen, saturation is largely responsible for the weak 
dependence of Ra(a2 on x p 

and,. ^Ba,^^ depends of 
R2i A2 can be described by the ratio between the inclusive nuclear structure functions was 
suggested for the case of an A dependent Pomeron structure and A independent flux. We 
have checked this conjecture using the results from Ref. jsjj and have found that it fails, 
the inclusive ratio being larger than the diffractive one. 



IV. SUMMARY 



Diffractive physics in nuclear DIS experiments could be studied at the electron-ion collider 
eRHIC. Hence it is interesting to extend the current ep predictions to the corresponding 
energy and targets which will be available in this collider. In this work we address nuclear 
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diffractive DIS and compute observable quantities like Udiff/<^tat and F®^ in the dipole 
picture. In particular, we have investigated the potential of eA collisions as a tool for 
revealing the details of the saturation regime. Since Udiff is proportional to a^ ip , diffractive 
processes are expected to be particularly sensitive to saturation effects. Moreover, due to the 
highly non-trivial A dependence of cr^p, diffraction off nuclear targets is even more sensitive 
to non-linear effects. Using well established definitions of (Tdiff and and a recent and 

successful parametrization of a dip-, we have studied observables which may serve as signatures 
of the Color Glass Condensate. Without adjusting any parameter, we have found that the 
ratio CT di ff/a to t is a very flat function of the center-of-mass energy W, in good agreement 
with existing HERA data. Extending the calculation to nuclear targets, we have shown that 
this ratio remains flat and increases with the atomic number. At larger nuclei we predict 
that approximately 37 % of the events observed at eRHIC should be diffractive. Moreover, 
we have analyzed the behavior of the diffractive structure function F®^ and found out that 
in certain regions of the (3 - xp space, the diffractive structure function F®^ decreases 
up to an order of magnitude when going from the lightest to heaviest targets. Finally, we 
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have found that, for nuclear targets, the contribution of the qqg Fock state becomes less 
important. 



Considering the results obtained in this paper and those presented in Ref. 24], we can 
conclude that eA collisions are very promising for the experimental confirmation of the 
non-linear effects of QCD. 
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